Dissipationless transport in low density bilayer systems 
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In a bilayer electronic system the layer index may be viewed as the z-component of an isospin-i. An 
XY isospin-ordered ferromagnetic phase was observed in quantum Hall systems and is predicted to 
exist at zero magnetic field at low density. This phase is a superfiuid for opposite currents in the two 
layers. At B = the system is gapless but superfiuidity is not destroyed by weak disorder. In the 
quantum Hall case, weak disorder generates a random gauge field which probably does not destroy 
superfiuidity. Experimental signatures include Coulomb drag and collective mode measurements. 
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In quantum well structures containing two separate 
two-dimensional electron gases in close proximity, an 
electron is described in terms of its position in the plane, 
its spin and its layer index. The latter can be regarded 
as an isospin-i, denoted by m, with the two layers being 
the two eigenstates of m^. States with spontaneous XY 
isospin-ferromagnetic order have been observed in quan- 
tum Hall systems at total Landau level filling factor 
1/ — 1 and were predicted to exist at B = for sufficiently 
low electron density. j2j The origin of isospin ferromag- 
netism is a favorable Coulomb exchange energy just as 
in the ordinary Stoner instability. 

For layer separation d — the isospin polarized phase 
breaks an SU(2) symmetry, and the problem maps onto 
the Stoner instability. [Q For d > 0, and in the absence 
of tunneling between the layers, there is an easy plane 
anisotropy since the direct Coulomb energy favors polar- 
ization in the XY plane ((toz) = 0) in order to avoid 
the cost of charge imbalance between the layers that oc- 
curs for {ttiz) 7^ 0. The angle of the magnetization m(r) 
relative to the x-axis is then described by a field if{r). 
Because the 'charge' conjugate to the phase ip is m^, the 
Goldstone mode |y] associated with the broken U(l) sym- 
metry at finite d corresponds to superfluid currents which 
are opposite in each layer. [^[j^J^ 

In this paper we study transport properties of the easy 
plane isospin ferromagnet, focusing on the effects of dis- 
order. At i? = we find that disorder weakens, but does 
not destroy the 'gapless isospin' superfluidity. The lack 
of time-reversal symmetry in the QHE case causes disor- 
der to induce a random gauge field which frustrates the 
system, but the Kosterlitz-Thouless transition probably 
survives weak disorder. The effect of random interlayer 
tunneling is a separate and different question. Q 

For _B = it is difRcult to quantify the range of param- 
eters, particularly rg, in which the isospin ferromagnetic 
state is the prevailing phase. It should lie between the 
low Tg paramagnetic range and the very high Vg range, 



where the system forms a bilayer Wigner crystal. For 
a single layer, Monte Carlo calculations find a fer- 
romagnetic transition in the rg ^ 20 — 30 range, and 
Wigner crystallization at « 37 ± 5 Q . The energy dif- 
ferences among the various possible phases are, however, 
very small so that a definitive statement is not possible 
- in fact, earlier calculations did not find a ferromag- 
netic transition. For a double layer system, Hartree-Fock 
(HF) theory predicts the existence of a broken symmetry 
isospin ferromagnetic phase. [|j HF tends to overestimate 
the stability of broken symmetry states, but its predic- 
tions are often qualitatively correct and such states fre- 
quently do occur at values of rg larger than predicted. 
While quantum Monte Carlo calculations are needed to 
obtain the precise density at which a _B = bilayer sys- 
tem will undergo the spontaneous isospin ferromagnetic 
transition, it is reasonable to assume, based on existing 
HF analysis Q|, that such a transition should occur at 
Ts w 20 — 30, a regime now realizable in hole systems. |^ 
An HF analysis of the isospin polarized phase starts 
with an Hubbard-Stratanovich decomposition of the 
Coulomb interaction, leading to the action 

[ dtdr\ij+idttp + —^Ij^HV - AaS']^^/j(r) 
J I 2m 

-p,(r)T/«(r)-p^(r)4^)(r) 

+ iv;,m2 + dr'n{r)Vg{r - r')n(r') 

-f iy"dr'm,(r)K(r-r')m.(r')} (1) 

In Eq. (|l|) ■0'*', Ip are fermionic fields for the electrons. The 
symmetric and antisymmetric densities are pg = ipa^o- 
and = S^^,i}}'^, with being the Pauli matrix for 
the ith-component of the isospin (i = x, y, z). The fields 
n(r), m(r) are auxiliary Hubbard-Stratanovich fields de- 
scribing symmetric and antisymmetric densities. We are 
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interested in the response of the system to a weak anti- 
symmetric vector potential A^, which is thus included in 
the action (a factor of | is absorbed in A^). In momen- 
tum representation, ^^(q) — ^ai'i) = 'iire'^d 
(for small q). For simplicity, we assume here that the 
true electron spin is fully aligned due to the Stoner in- 
stability and can be ignored (see however pO|]). 

In momentum representation, for small q, the sym- 
metric Hartree potential is vj^^ — Vs{q)n{q) while the 
antisymmetric is — Va{q)mz (q). The Fock potential 
Vex is approximated in Eq. (|l|) to be local, thereby ne- 
glecting the exchange contribution to the gradient terms 
which contribute to the isospin stiffness. We comment on 
the actual value of Vex and on consequences of its non- 
zero range below. In the system's response to A^ the 
symmetric field n(r) does not play any role and we omit 
it from following expressions. 

For fixed values of m and n, the action (|l|) describes 
non-interacting electrons under the influence of a space 
and time dependent scalar potential vjj \ vector poten- 

tial Aq and Zeeman field z + VexT^- In an x — y 
ordered state, the saddle point for the bosonic fields is 
n{r) — m2{r) — 0, and |m(r)| ~ A4, a non-zero constant. 
Conventional approximation schemes (HF, random phase 
approximation (RPA)) do not reliably obtain A4. Here 
we first assume full polarization {A4 — n, as predicted by 
HF), and later discuss the case of partial polarization. 

Due to an assumed lack of interlayer tunneling, the 
action (|l|) possesses a C/(l) symmetry. [|j Thus, in equi- 
librium the system picks an arbitrary direction for m. We 
write m = (cos((/7(r))x + sm{(p{r))y) + m^z, where (f 
is the angle between the planar component of the magne- 
tization and the ai-axis, and expect ip to be constant in 
the ground state and slowly varying in low energy exci- 
tations. The energy cost of a deviation from the equilib- 
rium magnetization is then expressed in terms of ip and 
niz, and should vanish for a uniform shift in <p(r). 

We now integrate over the fermionic fields and expand 
the familiar trlog{} term to second order in ip and m^. 
Within RPA the expansion is given in terms of the re- 
sponse functions Xi = ^{p^P^) and Xo = — (p^P^)- The 
effect of a slowly varying Aq on the x's can be separated 
out by means of a Gorkov approximation, where Aq is 
approximated not to vary in the range of r — r' and t — t' 
in which the response function is appreciable pl] |. The 
effect of A(j is then incorporated by the "minimal cou- 
pling" prescription iVtp — > (iV(p — Aq), and the response 
functions are calculated for Aa — 0. The RPA action is 
then (omitting the zeroth order term), 

Srpa^^ J duj y"dr|/9s|iVv3 - Aap 
+ —ml + 2xoVexM{Vex + 27re2d)m^(y9| (2) 

where 



Ps = - lim g-2(l + XyVexWexM^ 

q,cj— »0 
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- = - lim [1 + x.(2^e2(i + VexW^e'd + Vex) (4) 

The response functions XotXz,Xv are response func- 
tions of non-interacting electrons in a Zeeman field 
VexMx. For small q_,ijJ, 
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where A = MVex is the energy cost for flipping a spin, 
and the value of T) is discussed below. The U{1) invari- 
ance of the problem is the reason for 1 -I- Xy^ex being 
0{q^). The equation of motion for ruz, derived from (H), 
is the Josephson type relation (p = 2i:e'^dmz- 

The integral over rriz can now be carried out, resulting 
in an action in terms of Lp and A only, which is more 
transparent in space and time representation: 
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Eq. (^) is the action of a two-dimensional superfluid, with 
MT) being the superfluid "spin stiffness". If AiV > 0, 
the bilayer system responds to the vector potential A^ as 
a superfluid. and an antisymmetric current flows without 
dissipation. Eq. (^ reveals the existence of a longitudi- 
nal Goldstone mode that carries antisymmetric density 
and satisfies the dispersion relation 



Within RPA, the response functions Xy^Xz are 
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-} (8) 
in J 



Ilu + A + ea - ef3 + ir] lu - A + ea - ep + irj 

where |q:), are single particle eigenstates of the spin- 
independent non- interacting Hamiltonian, eQ,,e^ are the 
corresponding single particle energies, Pq is the density 
operator and /(e) is the Fermi function. 

Setting Lu — and expanding to second order in q, we 
find that for a clean system with full isospin polarization 
(2 A > /i, /i being the chemical potential): 
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This energy cost is the sum of single particle energies of 
eigenstates |k) of electrons in a Zeeman magnetic field 
that precesses in space in a constant rate Vtp and is 
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composed of one part originating from the anti- 

symmetric current induced by the precession of the field, 
and a second part, — ^^'^^^ , which reflects the slowing 
down of the symmetric motion due to the field preces- 
sion. There is no galilean invariance for antisymmetric 
currents so AdV ^ n/m. 

The disorder potential can be separated into symmet- 
ric and antisymmetric parts. The symmetric part affects 
2? much like non-magnetic disorder does in a conventional 
superconductor. For weak symmetric disorder (kpl 3> 1, 
and hence At ^ 1), the disorder- averaged matrix ele- 
ments in (ph are, 



|(«|Pq|/3)|' = -^ 



(10) 



where e = \{^a + c/s), D is the diffusion constant, v is 
the density of states, and |ea ~ c/jI < 7- Substituting in 
and paying attention to the dependences of y and D 
on I we find that the effect of symmetric disorder on the 
spin stiffness (|^) is of order l/Ar. 

Antisymmetric disorder modifies the capacitive energy 
term in (|l|) to be 2TTe^d J dr[mz(v) — mz^dis{T^)]'^ with ran- 
dom mz,dis- As is known from studies of, e.g., Josephson 
junction arrays, such a randomization in the equilibrium 
distribution of ruz reduces the superfluid density and can, 
if strong enough, induce vortex-antivortex pairs destroy- 
ing the superfluidity even at zero temperature. Here, 
since there are gapless Fermi surface excitations even in 
the superfluid, the resultant disordered phase may pos- 
sibly be a normal Fermi liquid with no long-range inter- 
layer phase coherence. 

Realistically, the disorder potential is made of com- 
parable symmetric and antisymmetric components. For 
weak disorder, then, antisymmetric currents flow with- 
out dissipation, although the superfluid density is sup- 
pressed. As the disorder gets strong, it eventually de- 
stroys the superfluidity. 

The superfluid spin stiffness A^P is also suppressed 
by finite temperature. Just as in an ordinary super- 
conductor, its temperature dependence originates both 
from the Fermi functions in (^) and from thermal fluc- 
tuations of vortex-antivortex pairs in (^(r). The spin 
stiffness, and with it long range order and antisymmet- 
ric dissipationless transport, disppear entirely above a 
Kosterlitz-Thouless (KT) transition temperature, whose 
precise value depends on both effects. 

An experimental probe of superfluidity of antisymmet- 
ric currents is the transresistance, or drag resistance, de- 
noted by Pd- In a drag measurement a current Ii is 
driven in one of the layers, while no current is allowed 
to pass through the second layer {I2 — 0) which devel- 
ops a voltage V2. Then, = — V2//1. For two identical 
layers, is the difference between the symmetric and an- 
tisymmetric resistances. In our case the latter vanishes. 
Thus, the transresistance equals the symmetric one, and 



the voltages on the two layers should he equal in magni- 
tude and direction. Since the superfluidity disappears at 
the K-T transition temperature, po would go down with 
increasing temperature. Note that for weakly coupled 
Fermi liquid bilayer systems is opposite in sign to the 
intra-layer resistance, and its magnitude increases with 
temperature. If the superfluid mode is lost due to disor- 
der, the antisymmetric resistance becomes appreciable, 
and the sign of po presumably becomes opposite to that 
of the intra-layer resistance. 

The excitation of the sound mode is another exper- 
imental probe. In the absence of isospin ferromagnetism, 
a double layer system has an antisymmetric acoustic plas- 
mon mode, which is overdamped by disorder as g ^ 0. 
|l2| Here, however, the sound mode (0) is an under- 
damped Goldstone mode. A density sweep experiment 
through the transition will therefore exhibit a sharp mode 
at low density which will get overdamped (at long wave- 
lengths) above the transition density. Another distinc- 
tion between these two collective modes is their behavior 
when d ^ Q. In that limit the Goldstone mode will have a 
long wavelength quadratic q^ dispersion, since the U{1) 
isospin symmetry changes into an SU{2), whereas the 
normal acoustic plasmon mode tends toward the single 
particle dispersion vpq- 

So far we have taken the exchange Fock potential to 
be local, and employed RPA. The U{1) symmetry of the 
approximate actions (^) and (^ is exact. However, other 
features of our analysis are not, of which we expect two 
to be most important. First, as a consequence of its finite 
range, the Fock potential renormalizes the dispersion re- 
lation of the electrons e(fc). The sum (|^) should then be 
evaluated with the renormalized energy dispersion, lead- 
ing to the replacements ^ — > Jp- and ^ — > ff in the 
intergal in (^, and affecting the spin stiffness. 

Second, it is conceivable that the Stoner phase is only 
partially isospin polarized, in contrast to the HF predic- 
tion of full isospin polarization. Interestingly, for a par- 
tially polarized state, and in the absence of disorder and 
electron-electron interaction, the g^-term in the sum (H) 
vanishes (i.e., V = 0), due the constant density of states 
h' characteristic of two dimensions. Spin stiffness is then 
induced by the deviation of v from a constant, caused 
by the renormalization of the energy dispersion by inter- 
action. Similarly, in the presence of symmetric disorder, 
V oc [{Di')' {fj.+A/2) — (Di/)' {fi~A/2)], where a prime de- 
notes differentiation with respect to energy. Again, the 
energy dispersion must deviate from parabolic for V to 
be non-zero. 

The physics of the isospin ferromagnet at filling factor 
1/ = 1 in the QHE regime is quite different from that 
at i? = 0. For two uncorrelated layers each at filling 
factor I' = 1/2, there will be gapless 'composite fermion' 
excitations In the presence of interlayer phase coherence 
however, the finite isospin stiffness leads to an energy gap 
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for symmetric excitations and a QHE plateau. 

Because of the energy gap for symmetric excitations, 
the fermions can be rehably integrated out to yield 

an Euclidean action which is a functional of the isospin 
unit vector a{r) = m(r)/n(r). The Hamiltonian density 
corresponding to that action is 



+ V,{r)6n{r) + + <5n(r)]a" (r). 



(11) 



r is the double layer capacitance per unit area (including 
Hartree and exchange contributions), Vs^a are the sym- 
metric and antisymmetric parts of the disorder potential 
and jO ~ IK is the exchange-induced spin stiffness. |^,|| 

A crucial consequence of the quantized Hall conductiv- 
ity is a constraint relating the symmetric fermion density 
n(r) = no + 6n{r) to the topological (Pontryagin) density 
of the isospin field a: \ 



6n{r) = ( -^a^y ) -LgM^e^,^ ^-^^^^3,^^ (12) 



where Roman indices range over x,y, z and Greek range 
over x,y, and no is the average density. 

Taking advantage of the easy-plane anisotropy and 
noting that the XY phase angle field ip contains vortex 
singularities, we integrate out the massive fluctuations 
and find that the lack of time-reversal symmetry causes 
the disorder potential to generate a gauge field yielding, 
in the high temperature classical limit, a 2D XY model 
with random Dzyaloshinskii-Moriya interaction 
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[XVaiTi,)~Vs{R,)Qj]M, (13) 



where Qj = ±1 is the vorticity of the jth vortex ('meron' 
[^,H), Mj = ±1 is a flavor index indicating the sign of cr^ 
in the vortex core, A is a non-universal constant related 
to the core size, and we have dropped various irrelevant 
terms (e.g., a random contribution to p). The gauge field 

^ = 4^1?^^+' where J+ = --^(Jxy'^'^'d^Vs is propor- 
tional to the symmetric Hall current. The TVa term is the 
local density imbalance. If there is a Hall current flowing 
when there is a density imbalance, then more current is 
flowing in one layer than the other and the superfluid 
mode |^,|| J_ ~ p(V(^ -I- a) is necessarily excited. This 
is the physical interpretation of the gauge potential a 
which causes these currents to flow. The phase field tp 
in ( |l3| ) contains both the singular and smooth parts, and 
thus the first term in (^|) mediates a logarithmic inter- 
action between the merons. Unlike the SU(2) symmetric 
case, where a symmetric disorder potential introduces a 
random gauge field [ p^ , in our case both symmetric and 
anti-symmetric parts are needed. 

The gauge field a is random with a finite correlation 
length. Ensemble averaging over a closed contour dj of 



perimeter L gives ( a • dr) = 0, and ( [ a • c?r] ) ^ 
LP with Q = \. This is the gauge glass model for which 
it is known that the KT transition is destroyed in the 
limit of strong disorder, jl^ For weak disorder the phase 
diagram has proven difficult to determine, [ p^ but it is 
likely that the KT transition survives. This result can be 
motivated by noting that in order to have an isolated flux 
quantum through a single plaquette, the vector potential 
would have to fall off like l/r and it follows that the 
random potential Vs.a would have to diverge. 

To conclude, our results call for an experimental search 
for the i? = phase coherent state in large Ts bilayer hole 
systems, and for a study of how superfluidity in this phase 
and its analogous phase at = 1 is affected by disorder. 
Superfluidity could be detected either in light scattering 
from the Goldstone mode or by drag transport. 
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